Interaction of electromagnetic perturbations with infalling observers inside spherical 

charged black holes 
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The electromagnetic radiation that falls into a Reissner-Nordstrom black hole is known to develop 
a "blue sheet", namely, an infinite concentration of energy density at the Cauchy horizon. The 
interaction of these divergent electromagnetic fields with infalling matter was recently analyzed (L. 
M. burko and A. Ori, Phys. Rev. Lett. 74, 1064 (1995)). Here, we give a more detailed description 
of that analysis: We consider classical electromagnetic fields (that were produced during the collapse 
and then backscattered into the black hole), and investigate the blue-sheet effects of these fields on 
infalling objects within two simplified models of a classical and a quantum absorber. These effects 
are found to be finite and even negligible for typical parameters of a supermassive black hole. 



I. INTRODUCTION 

The Kerr-Newman black hole (BH), which in view of 
the no-hair theorems is expected to be the stationary 
outcome of gravitational collapse, is an exact (electro-) 
vacuum solution of the Einstein field equations. However, 
one expects a generic collapse process to be accompanied 
by perturbations, which may exist before the onset of the 
collapse, or develop during it. Consequently, one would 
not expect to find an astrophysical BH to be an exact 
Kerr-Newman solution, but rather a perturbed one. 

The possibility to fall into a black hole and re-emerge in 
another universe is one of the most intriguing open ques- 
tions of General Relativity [jlj. The spacetime of unper- 
turbed BHs, such as Kerr, seems to allow this possibility. 
However, the Kerr geometry is highly symmetrical; It is 
not a priori clear whether more realistic solutions to the 
Einstein equations, which are not unperturbed, still al- 
low for the possibility to traverse the inner horizon safely 
and re-emerge in a different universe. Consequently, per- 
turbed BHs have been under investigation during the last 
three decades. 

Realistic BHs are the outcome of gravitational col- 
lapse. In such a realistic (a-symmetric) collapse, non- 
vanishing multipole moments of various fields develop in 
the star, and consequently electromagnetic and gravita- 
tional waves are emitted from the surface of the collaps- 
ing star. As these waves propagate outwards, some frac- 
tion of them is backscattered off the spacetime curvature 
and captured by the BH. This process leads to a "tail" 
of radiation, which at late times decays according to an 
inverse-power law both for the spherical case |^-|5|] and 
for the spinning case 

To re-emerge in another universe it is necessary to tra- 
verse a certain null hypersurface known as the inner hori- 
zon or the Cauchy horizon (CH). In order to cross the CH 
safely, it is necessary that the spacetime be either non- 
singular there or have at the most only a weak singularity. 
Therefore, it is troubling that the CH is a surface of in- 
finite blue-shift, i.e., infalling radiation (electromagnetic 



or gravitational), even very mild and well-behaved in the 
external universe, is infinitely blue-shifted at the CH M. 
There are two types of problems which may cause difficul- 
ties for an observer who wishes to cross the CH. First, a 
divergent spacetime curvature develops on the CH. This 
curvature singularity results from two sources: the infi- 
nite blue shift of infalling gravitational waves (which lead 
to the divergence of the gradients of the metric perturba- 
tions) and the infinite energy-momentum tensor (which 
can be taken as a dynamic source term for the Einstein 
equations) of the infinitely blue-shifted infalling electro- 
magnetic radiation. Second, the same infinitely blue- 
shifted electromagnetic waves cause an infinite flux of 
radiation, which might heat any infalling observer un- 
limitedly. 

Penrose || argued that the CH was unstable against 
small perturbations. His arguments for the infinite blue 
shift at the CH relied on a geometric-optics approxima- 
tion. Penrose argued that perturbations originated in an 
infinitely long (external) time are concentrated in a finite 
(proper) time of the infalling object. Thus, unless these 
perturbations decay at least exponentially fast in exter- 
nal time, the infinite concentration will lead to the blue 
sheet. Later, other works considered the wave equation 
for the evolving perturbations. Simpson and Penrose (Tc[ | 
re-affirmed numerically (for linear electromagnetic per- 
turbations) the qualitative arguments of Penrose. Gursel 
et al and Chandrasekhar and Hartle [Q calculated 
the projection of the energy-momentum tensor on the 
world-line of an infalling observer, and inferred from its 
divergence on the CH that the radiation absorbed by the 
observer was also divergent. Nevertheless, the fundamen- 
tal fields (i.e., the scalar field, the electromagnetic four- 
potential, and the metric perturbations) were found to 
be regular. It is the gradients of the fundamental fields 
which diverge on the CH. Therefore, it was suggested 
by Ori Jl3|,[l4| that although there is a true curvature 
singularity at the CH, this singularity is rather weak. 
Namely, despite the divergence of the tidal force, the ac- 
tual tidal distortion experienced by infalling observers 
(as they hit the singular CH) is finite, and for typical pa- 
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rameters - even negligible. To cross the CH safely there 
still remains, though, the other aforementioned potential 
problem. Namely, the possible complete burning up of 
any physical body at the CH, due to the divergent elec- 
tromagnetic field there (and the associated energy flux). 
This subject was considered recently by Burko and Ori 
fL5| . In this Paper we give a more detailed account of 
that work. 

Physical BHs are expected to spin very fast. However, 
it turns out that the mathematical analysis involved with 
the evolution of the multipole moments is very compli- 
cated, due to the axial symmetry of the Kerr background. 
For this reason, it is often common to work with a toy 
model, within which the mathematical analysis is much 
simpler. Of course, the toy model should preserve the 
most essential and relevant properties of the realistic BH. 
Thus, most work (including the present one) is done in 
the framework of the Reissner-Nordstrom (RN) BH. (The 
RN solution is involved, however, with a certain compli- 
cation, which arrises from the coupling of the gravita- 
tional and the electromagnetic fields.) The RN spacetime 
is the unique electrically charged, spherically symmetric, 
static vacuum solution of the coupled Einstein-Maxwell 
equations. For obvious reasons physical BHs are not ex- 
pected to be significantly charged. The vanishing angular 
momentum of the RN BH is also an unrealistic feature. 
Nevertheless, it turns out that the internal causal struc- 
ture of the RN BH is very similar to the internal struc- 
ture of the Kerr BH: In both spacetimes the singularities 
are timelike, and are located beyond a CH; Both space- 
times have a wormhole-like topology, which may allow 
for a travel to other asymptotically- flat universes. Con- 
sequently, it is believed that the RN BH is a physically 
justifiable model for realistic BHs. 

In this paper we analyze the interaction of an infalling 
object with the divergent electromagnetic field which we 
expect to find at the CH. Throughout this paper we shall 
assume that the infalling object is much smaller than the 
typical radius of curvature near the CH. However, elec- 
tromagnetic perturbations in the RN geometry are al- 
ways accompanied by gravitational perturbations. Typi- 
cally, the latter produce a diverging curvature at the CH 
(due to the infinite blue shift), which means that the ra- 
dius of curvature vanishes there. Our assumption is valid 
only if — for the sake of evaluating the radiative electro- 
magnetic effects — we ignore the gravitational perturba- 
tions. The modification of the radiative interaction by 
the metric perturbations is obviously a non-linear effect, 
as it is quadratic in the perturbation's amplitude. This 
non-linear effect remains the subject of future research. 
We do not expect, however, this non-linear effect to sig- 
nificantly alter the linear-order interaction. The gravita- 
tional analogue of this problem — the object's interaction 
with the divergent tidal forces — demonstrates this rea- 
soning: As implied from the analysis of Ref. O] on the 
strength of the CH singularity in spinning BHs — where it 



has been demonstrated that the non-linear gravitational 
interaction with an object may be negligible — we do not 
expect higher-order contributions to change the general 
picture significantly. Thus, in this work we restrict our- 
selves to linear effects only, and take the background to 
be RN (and not a gravitationally-perturbed background). 
In addition, recent fully-nonlinear numerical simulations 
have shown that at the asymptotic past of the CH the 
metric perturbations vanish, in accord with perturbation 
analyses p6| . 

The organization of this Paper is as follows: In section 
II we summarize the definitions and the notation we use. 
In section III we briefly review the formalism given by 
Chandrasekhar jl7) for the determination of the tetrad 
components of the electromagnetic field tensor F a p at the 
CH for given perturbations. In section IV we obtain the 
tetrad components of the electric and magnetic fields for 
initial perturbations which decay according to an inverse- 
power law. In section V we transform these tetrad com- 
ponents to their tensorial counterparts and write them 
in the rest-frame of a freely- falling observer. In section 
VI and section VII we use these fields to calculate their 
interaction with (very simplified) classical and quantum 
absorbers, respectively. We then discuss (section VIII) 
the results, and argue that although the electromagnetic 
fields diverge at the CH, the interaction of the field does 
not necessarily cause ultimate destruction of the infalling 
observer. 

It should be noted that we only treat here classical 
radiation, and ignore quantum processes, especially pair- 
production effects. When these effects are taken into 
consideration ]l8[ ], it should be expected that they may 
change our results here considerably. Yet, we believe that 
our main conclusion, namely, that the singular CH is not 
the edge of spacetime, will still be relevant even after 
consideration of the quantum effects. 

II. DEFINITIONS AND NOTATION 

We write the line element of an unperturbed RN BH 
in the form 

ds 2 = e 2v (dx ) 2 - e 2 ^ (dx 2 ) 2 -r 2 dn 2 , (1) 

where the coordinates are (x° x 1 x 2 x 3 ) = (t <p r 6), dQ 2 
is the unit two-sphere line element, and the metric coef- 
ficients are e 2v = e~ 2 ^ 2 = (r 2 - 2Mr + Ql)/r 2 = A/r 2 , 
M, Q* being the mass and electric charge, respectively, 
of the RN BH, and where r is the radial Schwarzschild 
coordinate, defined such that circles of radius r have cir- 
cumference 2irr. The general form of the line element (1) 
is preserved under polar perturbations (sometimes called 
even-parity perturbations); On the other hand, axial per- 
turbations (called also odd-parity perturbations), will 
lead in general to non- vanishing off-diagonal metric. [Ax- 
ial perturbations are characterized by the non-vanishing 
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of the metric functions cj,<72,93 (the non- vanishing of 
these metric coefficients induce a dragging of the inertial- 
frame and impart a rotation to the BH) , while polar per- 
turbations are those which alter the values of the metric 
functions v, fj,2, ^3 and ip (which are in general non-zero 
for the unperturbed BH).] Therefore, the form of the met- 
ric of a generally-perturbed RN BH will be more compli- 
cated than the line element (1). It has been shown Jl7| , 
that a metric of sufficient generality is of the form 



ds 2 



e 2lJ (dx ) 2 - e 2 <" ( dx 1 - cj dx° - q 2 dx 2 
q 3 dx 3 ) 2 - e 2 " 2 (dx 2 ) 2 - e 2 ^ (dx 3 ) 2 . 



(2) 



Since the unperturbed RN background is spherically 
symmetric, we may consider only axisymmetric pertur- 
bation modes without any loss of generality. (This is be- 
cause all non-axisymmetric modes can be received from 
the axisymmetric modes, if the unperturbed spacetime 
is spherically symmetric Jl7|.) The metric (2) involves 
seven functions, namely, v, ip, /X2, /X3, u>, q%, and q%. Bea- 
cuse the Einstein equations involve only six independent 
functions, not all of the seven functions can be deter- 
mined independently, and there is a gauge-fixing freedom 
on the metric coefficients. However, this gauge freedom 
involves only the metric coefficients u>,q2, and 53; There 
is no gauge freedom in the determination of the polar 
perturbations. 

The horizons of the RN BH are the event horizon r + 
and the inner horizon r_ , which are located at the roots 
of A, namely, at r± = M ± (M 2 - Ql) 1/2 . We de- 
note the surface gravity of the event horizon and the 
CH by k± = (r+ — r_ ) jr\ , respectively. We define the 
Eddington-Finkelstein null coordinates u = r* — t and 
V — 7'* + t, where is the Regge- Wheeler 'tortoise' coor- 
dinate defined by dj dr* = (A/r 2 )d/ dr. The coordinate 
t is spacelike between the event and the Cauchy horizons, 
and we take t = +00 at the event horizon. In this Paper 
we are interested in the sections of the event horizon rep- 
resented by u — — 00 and of the inner horizon represented 
by v = +00. (These are the sections which intersect in 
the standard Penrose diagram at future timelike infin- 
ity of the external universe.) We assume that the object 
moves along a typical radial world line that intersects the 
event horizon and the CH at some finite values v = vq 
and u = uq, respectively. Accordingly, the trajectory of 
the object can be described by the function r(r) and by 
Mo, where r is the proper time of the infalling object. We 
set r(r = r_) = 0. 



III. THE CHANDRASEKHAR FORMALISM 

In this section we briefly review the Chandrasekhar 
formalism for the evolution of the polar perturbations 
p7[ , and describe the algorithm which can be constructed 
from it for the determination of the Maxwell tensor. 



The Chandrasekhar formalism fails for the description 
of dipole perturbations, i.e., for the I = 1 modes in the 
multipole expansion of the perturbations. However, for 
the treatment of all modes with multipole order higher 
than the dipole, it can still be used. The formalism for 
dipole polar modes was treated by Burko and is given 
in Ref . jl9| . The modified formalism of Ref. can be 
used in an analogous way for the determination of the 
dipole perturbations. 

Let us consider a RN BH perturbed by the electromag- 
netic waves, which we expect to exist for any charged BH 
created by a collapse process. Let the metric of the un- 
perturbed RN background be given by Eq. (0). The 
metric form (0) does not change when the collapse is 
endowed with polar perturbations jl7). [Axial pertur- 
bations, on the other hand, will lead to the appearance 
of other (off diagonal) non- vanishing metric coefficients, 
and thus complicate our calculations considerably. For 
this reason, and as knowledge of the perturbations of 
either of the classes enables us to know the perturba- 
tions of the other class (to be explained below) , we con- 
sider here polar perturbations only. It should be remem- 
bered that our entire analysis is linear in the perturba- 
tions, and therefore effects of non-vanishing off-diagonal 
metric coefficients are negligible.] We take — after Price 
p|, who considered the uncharged Schwarzschild back- 
ground, and Bicak Q, who extended the treatment to 
the RN background — the perturbing fields to decay ac- 
cording to the (21 + 2) inverse-power in external time. 
The infalling radiation propagates through the curved 
spacetime of the BH. Due to the RN background, the 
(linear) gravitational and electromagnetic perturbations 
are coupled. However, their equations can be decoupled 

The metric perturbations, analyzed into their normal 
modes with a time-dependence e lCT * are governed by the 
Moncrief-Zerilli [^TJ equations, which are a pair of decou- 
pled one-dimensional wave equations of the form 



2y(±) 



dPZ. 



dr 2 



a 2 Z, 



(+) 



(3) 



Here T^ (±) = ±ft dfc/ dr*+/3 2 ff+nfr, where k = /i 2 (/J 2 + 
2), fa = qj and f t = A/[r 3 (/j 2 r + q/)}. Here, ij = 1,2 
and i ^ j. We also take fi 2 — (I — I) (I + 2), <ji, 2 = 
3M± a/9M 2 + 4Q 2 /?, I being the multipole order of the 
perturbing radiation. The superscript (+) denotes polar 
perturbations, and (— ) denotes axial perturbations. As 
shown in Ref. knowledge of modes of each of the 
parities enables us to find the modes of the other parity, 
so without any expected loss of generality we shall restrict 
ourselves here to the treatment of polar perturbations 
only. (It should be stressed that these expressions are 
inapplicable for dipole modes. For the determination of 
the latter one must use the formalism of Ref. 

Through the end of this section we shall briefly sum- 
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marize the algorithm we can construct from the Chan- 
drasekhar formalism (or from its generalization to dipole 
modes) for the calculation of the perturbed metric func- 
tions and the components of the Maxwell field strength 



tensor, given the initial perturbing fields Z 



(±) 



1,2). 



When we give refernce to a specific equation in Chap- 
ter 5 of Rcf . jl7| , a letter "C" will precede the equation 
number. We first find the functions by the algebraic 
equations (C186)-(C187). Now, we calculate the function 
$byEq. (C196). The next step is to obtain the functions 
defined by Eqs. (C190), and complete the solution with 
Eqs. (C191)-(C195). The perturbations of the metric (1) 
will be given then by Eqs. (C158) and the perturbations 
in the tetrad components of the Maxwell field strength 
tensor will be given by Eqs. (C159-C160). In the follow- 
ing we shall transform from the tetrad components of the 
Maxwell tensor to their tensorial counterparts using the 
tetrad 

e m = e v {\ 0) 
e (1)Al = rsin(9(0 1 0) 
e (2 )„ = e- v (0 1 0) 
e (3)M = r(0 1). 



IV. DERIVATION OF THE MAXWELL FIELD 
NEAR THE CH - TETRAD COMPONENTS 



The Components of B M „ 

Our goal now is to find an approximate expression [to 
the leading order in (r — r_) and (k^uq)^ 1 ] for the elec- 
tromagnetic field an infalling observer measures on ar- 
rival at the CH. We shall find that some components 
of the Maxwell field strength tensor diverge there, while 
other components remain finite. Therefore, we shall re- 
strict ourselves here to the evaluation of the divergent 
components only. (The other components can be found 
analogously.) We shall consider here an I > 2 polar mode 
of infalling radiation. It should be stated that the most 
dominant effect is not caused by the / > 2 modes but by 
the / = 1 mode. The reason we choose here to discuss 
the I > 2 modes rather than the I = 1 mode is that the 
formalism of Ref. JItJ is inapplicable for the treatment 
of dipole modes, as mentioned in section III (see Ref. 
|l9| for details). It turns out, that when one calculates 
the dipole perturbations according to the generalized for- 
malism of Ref. Jl9) the results remain qualitatively un- 
changed. Therefore, we may consider here only the I > 2 
modes. When we conclude the perturbation analysis we 
shall give the final result for the dipole perturbations too. 

Let us take, then, the electromagnetic and gravita- 
tional perturbations near the CH to be 1 1 1 : 



-(21+2) +bu -(2l+2) 



(4) 



In what follows, we consider an isolated charged BH, 
surrounded by electromagnetic waves, which we treat 
as a linear perturbation. (In fact, because of the non- 
vanishing electric field of the background, this linear per- 
turbation consists of both electromagnetic and gravita- 
tional waves pl]Jl2]] .) We calculate the asymptotic be- 
havior of the electromagnetic perturbation near the CH. 
The class of perturbations that we consider here is the 
one which is inherent to any non-spherical gravitational- 
collapse; These are the electromagnetic perturbations 
which result from the evolution of non-vanishing elec- 
tromagnetic multipole moments (in the star) during the 
collapse. When these perturbations propagate outwards, 
some fraction of them is backscattered off the spacetime 
curvature and captured by the BH. This process leads to 
a "tail" of infalling radiation at the event horizon which 
at late times (v ^> M) decays like (v/M)~^ 2l+2 \ where I 
is the multipole order of the mode [g| . We treat this elec- 
tromagnetic field according to the formalism by Chan- 
drasekhar JItJ for I > 2 polar modes and the extension 
of this formalism by Burko jl9| for dipole (1 = 1) polar 
modes. (The extension of the formalism is needed as we 
a posteriori find that the effects we study are dominated 
by the dipole mode. This mode is not treated properly 
by Ref. 0.) 



(5) 



We now use the algorithm described in section III to 
obtain the electromagnetic field near the CH. From z[ + 2 

we can find the functions by Eqs. (C186-C187). 

The function $ is given by Eq. (C196). Substituting Eq. 
(C196) in Eqs. (C186-C187), we obtain 



$(r,t) 



nr\qiq 2 \ 



'2nQ*qi 



/(+) 



}6 
dr. 
mr 



As we are interested in the development of the perturba- 
tions near the CH, we write an approximate expression 
for <3>(r, t) to the leading term in K_r„. To do this we 
expand in r — r_ , and obtain 



■ dr : 



-r„/2 



nr 2 _ + 3Mr- - 2Ql 



0{e 



') 



dr* 



Substituting in $ we get 
$(r*,f) k, a x 



e »' 



(a 2 a + a 3 c)v-^ 



(a 2 b + a 3 d)u-( 2l+ V 



dr* 



(6) 



where 
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a>i = {<& + \QiQ2\) 1 ( 
a-i = nr\f( 



3Mr- - 2Ql) 



\/{2njQ*qi 
a 3 = nrqi - \J (2n)yJ (-qiq 2 )Q*. 

It can be verified (for a formal proof see Ref. (||]) that 
the integral J°° e~ 9Z z~( 2i+2 ) dz can be represented by 
the asymptotic series 



9 



-gx x ~(2l+2) 



p=0 



. (P+2Z + 1)! 
(2Z + l)! 5 f 



Since we are interested primarily in very large values of x, 
i.e., we are interested in the regions of spacetime closest 
to the CH, the zeroth order approximation of the asymp- 
totic series suffices for our needs. Therefore, we shall take 

f°° e -^ z -(2'+2) dz ~ X -gx -{-il+2) Tq uge th j for th 

evaluation of $ we first change the variables in Eq. (|6|), 
and integrate the two terms separately. Thus, we obtain: 



*(r.,t)« -^e 



( a2 a + a 3 c)«- (2i+2) 
+ ( a2 & + a 3 d)u^ 2i+2 )] . (7) 

In order to calculate Bq3, we first need to have H^y. 
[See Eqs. (C190) and (C195).] We readily find that 

(+) _ dr* d (+ ) 



-(21 + 2) 



-2v 



moa\ 



(W 



-(21+3) 



f3u 



-(21+3) 



where a = q\a — (— q\q2) 1 ^ 2 c and (3 = q\b — {—qiq2) 1 ^ 2 d. 
We can obtain the following approximate expression to 
the dominant term in e~ v (as e~ v diverges on the CH it 
is clear that the stronger the dependence of the exponent 
in v the faster the divergence): 



#03 (r, t) 



(21 + 2)Q^i 
A (g 2 + 191921) 



-(2i+3) 



j3u 



-(21+3) 



(8) 



[In Eq. (||) we kept only the leading term in e~ v '.] Simi- 
larly, we also find the dominant term in e~ v of B23(r,t) 
near the CH. Namely, in view of Eq. (C190) and Eq. 
(C194) we get 



1 



9 1 

r q{ 



or 



-(21+2) 



+ f3u 



-(21+2) 



We thus see that both B 03 (r,t) and -B 23 (r, i) are the re- 
sults of linear differential operators acting on the per- 
turbing fields z\ + \ As shown in section III [see Eqs. 



(C159)-(C160)] the formalism can be now used to obtain 
frequency-dependent tetrad components of the Maxwell 
tensor. As the expressions for the functions B^ v are in- 
dependent of the frequency, it is clear that there is a need 
to adapt the two representations of the fields (the tem- 
poral representation and the frequency representation). 
It turns out that the tetrad components of the Maxwell 
tensor which lead to divergencies are i*(o)(3)(r, a) and 
F (2)(3)(r, a). 

The Components of F(jA( u \ (r, t, ff) 

Our goal now is to obtain the frequency-independent 
expression for the electromagnetic field. As -F(o)(3) ( r , &, 9) 
does not involve the frequency [in view of Eq. (C159)], 
it is possible to compute -F(o)(3) {r, t, 9) directly, without 
considering the subtleties of the Fourier transform. It 
can be shown that for the correct performance of the 
transformation from the frequency representation to the 
temporal representation one should only replace the fre- 
quency a of Eq. (C160) with —id/dt. (For a rigorous 
proof of this scheme — which is not as trivial as it may 
seem — see Ref. p2fl-) We thus obtain: 



F (om (r,t,9)^ - 



(l + x)(j, e- 
9? r " 



31 92 1 r 



-(2Z+3) 



(9) 



F (2)(3) (r,i,0)«(Z + l) 



lie 



1 



r q\ 



-(21+Z) 



(3u 



-(21+3) 



I Si 92 1 
P-2.B. 



(10) 



Here, P/(cos#) denotes the Legendre polynomial of order 
I. 



V. DERIVATION OF THE MAXWELL FIELD 
NEAR THE CH - TENSORIAL COMPONENTS 

This section is built in the following way: We first 
write down the tensorial components of the Maxwell ten- 
sor in the Schwarzschild coordinates, which are known 
(see, e.g., Ref. ^3|) to be singular at the CH. We then 
transform from the Schwarzschild coordinates to Kruskal- 
Szekeres coordinates, which are regular at the CH. Fi- 
nally, we transform from the Kruskal-Szekeres coordi- 
nates to the rest frame of the infalling object. (The mo- 
tivation behind this last transformation is given below.) 



The Schwarzschild Coordinates 

After finding the tetrad components of the electro- 
magnetic field we shall now find the tensorial compo- 
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nents. From now on we omit the explicit angular de- 
pendence of the Legendre polynomials, to make the ex- 
pressions simpler and more compact. We may do this 
because none of the transformations to be performed be- 
low changes the angular coordinates. To allow for this 
omission we re-define the components of the Maxwell 
tensor. Schematically, we separate the variables by 
F( M )(u)(r,t,9) = F^)^(r,t)Q^ v )(9). Hence, from now 
on Fr^)r v -j(r,t) should be understood accordingly. The 
tensorial components can be obtained from the tetrad 
components by the transformation = F^^e^ e[f\ 
namely, 



F 03 (r,t) 



-re u F m3) 
9? + 1 9192 1 



(r,t) 



(11) 



Now, we define the Kruskal-Szekeres future directed null 
coordinates U, V by: 



In 



In 



1 



K-V 



1 



:K-U = 



(15) 



(16) 



To obtain the form of the metric in these coordinates we 
first need to have an explicit expression for (r), because 
in transforming from the the Eddington-Finkelstcin co- 
ordinates to the Kruskal-Szekeres coordinates we find for 
the line element 



ds 2 



{r+ r){ { r -h K - r > M dU dV - r 2 (U, V) dQ 2 . 



F 23 (r,t) = re-»F i2){3) (r,t) 
e- 2v (l + l)n 



9? 



I 9i 92 | 



-{21+3) 



f3u 



-{21+3) 



(12) 



Transforming to Regular Coordinates 

Transforming to Kruskal-Szekeres Coordinates 

Now, we wish to transform the Maxwell strength field 
tensor to a co- moving reference frame, i.e., a frame in 
which the infalling observer is at rest. We do that in order 
that we could relate the results for different observers and 
because of the fact that to use a (classical or quantum) 
local theory for the matter-radiation interaction we need 
to express the electromagnetic field as measured by the 
physical system in question, and as a function of its local 
(proper) time. First, we transform from the coordinates 
(t (/> r 9) to the coordinates (t <fi r* 9). It is clear, that the 
only component (out of the two relevant components) of 
F^y which is changed by this transformation is F 23 . In 
fact, we get that 



F r 
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Now, we transform to the coordinates (u <f> v 9). 
(The coordinates u and v arc defined in section 
II.) The line element for the unperturbed RN back- 
ground in theEddington-Finkelstein coordinates is ds 2 = 
— ^rdudv — r 2 (u, v) dQ 2 . We thus find that 



F U — F r — F t g — — -tj 

9? 



f±^t-Ru-^ (13) 



9i 92 1 



F vg = F r , 6 + F tg = f + ^. av-VW. (14) 
9J + |9i92| 



As we defined r* in section II only through its differen- 
tial, we realize that when we integrate to obtain r*(r) 
we may add an arbitrary integration constant. This in- 
tegration constant can be chosen in such a way, that the 
line element near the CH assumes a quasi-Minkowskian 
formQ, i.e., ds 2 = [-1 + 0(UV)]dU dV - r 2 dil 2 . Thus, 
we integrate dr*/ dr and obtairjj 

1 T-L. — V 1 T — T — 

7'* (t) = (r — r_ ) H In In 

k + r + k_ r_ 



1 



In- 



(r+ — r _)i+K-/K 4 
Inserting this in the line element, we find 

l+K— / K+ 



ds' 



r_ s • 
r 



r 2 (U,V) dil 2 . 



e «- ('-'■-) dUdV 
(17) 



Here, r is the implicit function of the coordinates UV — 
4exp[— K_r»(r([/, V))]/n 2 _. This metric is of course reg- 
ular in the domain between the two horizons, and at the 
CH, as should be expected from Ref . |^4|,^5| . Due to the 
regularity of the metric ([jl]) at the CH, we can work on 
a sufficiently small neighborhood where spacetime is as 
close to Minkowski spacetime as we wish (which is, in 
fact, trivial, as any non-singular spacetime possesses this 
property). In these coordinates we obtain: 



*This is the line element not only exactly on the CH but 
also very close to it. (Of course, this line element is regular 
at the CH.) We shall use this quasi-Minkowskian form when 
we construct interaction models for the radiation with matter 
(see sections VI and VII). 

* Notice the difference between this form of r*(r) and the 
forms of Ref. Also notice that in Ref. Ill the defini- 

tion for the surface gravity is different from ours. 
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(18) 



F UB = -^= 2 ^t 2) (19) 



K-U 



k-U 9i + |gig 2 | 



the proper time of the infalling object. It is clear, 
that these transformation equations define the special- 
relativistic Lorentz boost generating the co-moving coor- 
dinates uniquely. The solution of these equations is 



at ot 

8T dZ 

dz dz 

dT dZ 



T -Z 
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Transformation to the infalling object's rest-frame 



Eqs. ( |18| , |19[ ) are not enough for our needs. The reason 
for this is two-fold: First, we should like to calculate the 
interaction of the electromagnetic field near the CH with 
the matter comprising the infalling object. A description 
in a flat spacetime of course simplifies the intricate ra- 
diative processes considerably. A flat-spacetime descrip- 
tion also allows us to use the standard notions of elec- 
tric and magnetic fields. Second, The region where the 
electromagnetic field assumes exceptionally high values 
is a very 'narrow' region near the CH. In this narrow re- 
gion, because spacetime is approximately Minkowskian, 
spacetime curvature is negligible in the background. Ad- 
ditionally, it will turn out (see sections VI and VII) that 
we shall need to describe the interaction of the infalling 
object with the radiation field along the entire trajectory, 
and, in addition, we shall also wish to compare different 
observers arriving at totally different points on the CH. 
(The reasons for this will become clear in sections VI and 
VII.) To do that we cannot be satisfied with a completely 
local description of the CH, and therefore we transform 
now to co-moving coordinates, which allow us to com- 
pare different observers, as for all we set the proper time 
equal to zero on arrival to the CH. Let us now define 
the coordinates Z, T by: U = Z - T and V = Z + T. 
(Note, that Z is a timelike coordinate and T is a space- 
like coordinate.) The coordinates Z, T are Minkowski- 
like coordinates, yet not appropriate for the description 
of the object's rest frame, as generally we should find that 
the infalling object is in motion relative to this reference 
frame. 

We now transform to co-moving Minkowski coordi- 
nates z, t, adapted to the trajectory in question. Namely, 
we demand that at the intersection of the orbit with 
V = 0, the newly defined coordinates be such, that z = 
(and t = 1). In addition, we set z = t — at that inter- 
section point. This transformation is thus defined by the 
coupled equations 

(1 0) = x a ' = ^hrx? 



n 



a/3' 



r a' j3' dx a dx^ olQ 

y dx a dx$ y 



where r\ a ^ is the metric tensor of a Minkowski space- 
time and a dot denotes differentiation with respect to 



We now define the null coordinates u = z — i and v = 
z + t. (Note that both u and v vanish at the intersection 
of the trajectory with V — 0.) Transforming to these 
coordinates wc find 
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du Fu6 



It can be shown that 

du du 
~dV = ' ~dU 

Hence, we get that 

1 



V, 



dv 
~dV 



U, 



dv 
~dU 



Fve — Jj^ve, 

It can be readily shown ^ 
l/(2r)e- K - u °/ 2 and V(r-) 



FuB — ~^F\J6- 

that on the CH £/(r_) = 
= -2fe K -" o/2 . Using this, 
it can be shown that on the CH v — UV and u = 
VU — Ar/K—. After performing all the transformations, 
we end up with the required expression for the divergent 
component of the Maxwell tensor as measured by an in- 
falling observer in his rest frame. We substitute 



v = In 
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k-VU 91 + |9i92| 
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, (21) 
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where 



C' = 



C"' 



4(Z + 1)// 
Qi + ki92| 
4(l + l)/z 

9? + \qiQ2\ 



\2i+3 1 



The form in which we represented F u g might obscure the 
simplicity of its meaning. In fact, all we need is the value 
on the inner horizon, where u — 0. We readily find that 
on the CH 

rp 1 + 1 r -M a -(21+3) 

rue = — ; 2~n \P u o 

r qf + |?i92| 

From this expression the regularity of F a g on the CH is 
self evident. 



The final form for the electromagnetic field near the CH that 
an infalling observer measures 



typical observers r is neither vanishing nor divergent, and 
is typically of order unity. (On the other hand, In |k_t| 
diverges, and uq is typically very large too.) From now 
on, for the sake of brevity we shall call these electromag- 
netic components simply E and B, respectively. Thus, 
we conclude that 



C ( 1 
E = —B w In \k—t\ + —K-Uo 

K—T \ 2 



-(22+3) 



(22) 



It is clear from Eq. (|2^) that the electromagnetic field — 
as measured by the infalling observer — diverges. This 
conclusion agrees with the results of Refs. [g ; 
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The above analysis was done for I > 2 polar modes. 
The analysis can be repeated for I = 1 polar modes anal- 
ogously using the formalism of Ref. |L9]]. When this is 
done, it turns out that the electric and magnetic fields can 
still be expressed by Eq. (|2^). However, the correct ex- 

-8a( K _/2) 5 /(<Zir-). 



pression for C is now C'(l = 1) 



We may construe each point of the matter comprising 
the infalling observer as being located in the center of its 
own spatial coordinate system, i.e., we may take 2 = 
or, equivalently, v = u = i. Identifying the coordinate 
i with the observer's proper time r, we realize that v, u 
vanish on the CH. We get that F u g remains finite on the 
CH, while F^e diverges. 

We can write the required divergent components of the 
electric and magnetic fields near the CH as measured by 
the infalling observer in an orthonormal Cartesian tetrad 
frame in which the observer is at rest. As the tetrad 
introduced in section III is orthonormal, it is only natural 
to take here the same tetrad base. Thus, we get: 



F (2)fe) - F * ee h e (y) = F ^ e (v)^ 



= -Fge, 
r 



(t)(y) 



fj.6 _ 



-F, 
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As the tetrad base is Cartesian and orthonormal, it is 
easily shown that both ef^x and e*^ identically equal 
unityQ. We now denote Ey = Fig and B s = —F s g, 
where F^g w F E g « F^g. Here, x and y are or- 
thonormal tetrad components in the object's reference 
frame, directed in the dj d<f> and the dj 08 directions, 

— Bx ~ 



respectively. Hence, we find that Ey 



C 
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From these ex- 



'— (ln|K_r| + iK_u +ln|r|)" 
pressions for the divergent components of the electromag- 
netic field we may neglect the term dependent on f, as for 



*A two-dimensional flat Minkowskian spacetime is spanned 
by the vectors z and i. As these vectors are orthonormal, the 
claim is readily proved. 



VI. INTERACTION OF THE RADIATION WITH 
MATTER: CLASSICAL ABSORBER 

We now consider the interaction of the electromagnetic 
field (^||) with the matter comprising the infalling object. 
Here, we consider a classical object, and in the next sec- 
tion we consider a quantum object. In both the present 
and the consecutive sections we assume that the object 
is much smaller than the typical radius of curvature be- 
tween the event and the inner horizons, and hence the 
effects of curvature are negligible. This arrises from our 
assumption that non-linear effects will not have impor- 
tant effects on the radiative interaction of the field with 
the infalling object. (See the discussion in section I.) 
Consequently, we can imagine the object as being at rest 
in its locally co-moving Minkowski frame when an elec- 
tromagnetic impulse of the shape ( |22] ) comes from null 
infinity and interacts with it. In what follows we shall 
use the co-moving Cartesian coordinates defined in sec- 
tion V, but for simplicity we shall omit here the 'bars.' 
Namely, we shall use the coordinates (r x y z). 

Despite the flatness of spacetime in the observer's rest- 
frame, the interaction of the electromagnetic field with 
the matter from which the observer is made is extremely 
complicated. Therefore, we take a very simplified (toy-) 
model for the matter-field interaction, which still embod- 
ies the most essential properties of more realistic inter- 
actions. We take the matter to be composed of classical 
"atoms." Each "atom" is composed of two point-like op- 
positely charged particles, with charges ±e and masses 
/j,±, correspondingly. We denote the reduced mass by \i. 
(Do not confuse fi here, which is the reduced mass, with 
fi in the previous two sections.) The system interacts 
with an external force, which in our case is the Lorentz 
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force induced by the blue-shifted incoming electromag- 
netic field. This external force acts to change the separa- 
tion distance between the two particles of the "atom." We 
presume small deviations from equilibrium (to be justi- 
fied a posteriori) , and hence take a linear restoring force, 
i.e., F — —/iaj 2 X, where X is the deviation (of the parti- 
cles' separation) from equilibrium, and 10 is the resonance 
frequency. The dipole is chosen to be aligned in the dj dO 
direction (to allow for a maximum interaction with the 
field). We take the initial conditions to be X = and 
X = 0. The excitation of the system by the field may be 
characterized by X, X, and also by the total absorbed 
mechanical energy £ c . We shall show that all three vari- 
ables are finite, small, and for typical parameters even 
negligible. 

The equation of motion is 



fiX + [iu?X = eE(r), 



(23) 



where E(t) is the divergent component of the electric 
field (^2[). (The contribution of the magnetic field is ne- 
glected, as the ratio of the electric term to the magnetic 
term in the expression for the Lorentz force is propor- 
tional to the system's internal velocity X, which is taken 
to be small — a presumption which is justified a posteri- 
ori.) 

The solution of this equation (with our initial condi- 
tions) is 

X(t) = —-e-^ T [ E(T')e luJT ' dr' + c.c. (24) 

2lLU [I J _ T 

XM = --e~ WT [ E{T')e^ T ' dr' + c.c, (25) 
2 fj, J_ T 

where T is the time of infall from the event horizon to 
the CH. Calculating the sum of the kinetic and potential 
energies, we find for the total absorbed mechanical energy 



1 



1 



£ C ( T ) = -^u) z X A + -\xX 



E(r')e WT ' dr' 



-T 



(26) 



In all of the three expressions (|24j) , (E5J) , and (|26|) we need 
to calculate the integral 



/= / E( T ')e WT ' dr' . 



(27) 



The evaluation of this integral for typical parameters is 
done explicitly in Appendix A. The evaluation yields [Eq. 
(A27)], to the leading orders in (k_uq) _1 



C 1 1 
I(t = 0) » — — ~ ^ ( -k_«o 



(2Z + 2)k_ V 2 
C 1 . (\ 



-(21+2) 



-_r i \r- ua 
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(28) 



Substituting in the explicit expressions for X, X, and £ c 
we obtain, to the leading order in (k_uq) _1 



£o(o) 



1 C" ; 



2(2/ + 2) 2 k 2 _^ \nj \2 



/' 1 - I 1 77 K - M 



-2(21+2) 



21 + 2k_ \(i J \2 



2 k_uj \2 



, (29) 
(30) 
(31) 



We find that the strength of the excitation depends on uo. 
However, uq cannot be directly evaluated by an outside 
observer, who wishes to predict the excitation strength 
should he jump into the BH. Thus, it would be worth- 
while to express the excitation strength in terms of the 
external parameters of the problem. This can be done 
once we realize that duo/ dvo = — 1. (The proof is given 
in Ref. \22\.) We find that the infalling observer can 
increase |uo| by simply waiting outside the BH before 
jumping in and thus increasing the value of vq. For a 
sufficiently large |ito| we get that £ c (0), X(0) and X(0) 
are finite and arbitrarily small. 

We conclude that despite the divergence of the electro- 
magnetic field, the excitation of the classical system (and 
the energy absorbed) is finite, and becomes arbitrarily 
small, for late-time observers (large Vq). Therefore, the 
behavior of the charged classical system obtained by the 
above analysis is regular, and despite the divergence on 
the CH of the external force acting on the system, the 
energy absorbed by it is finite and negligible for a suffi- 
ciently large vq- 



VII. INTERACTION OF THE RADIATION WITH 
MATTER: QUANTUM ABSORBER 

As a quantum analogue to the preceding classical sys- 
tem we take (again) a very simplified model, which, we 
believe, captures the essential properties for the interac- 
tion we study and neglects all irrelevant details. Let us 
deal then with a non-degenerate quantum system obey- 
ing the Schrodinger equation. We shall evaluate the ex- 
citation of the system by considering the transition of 
the system from its ground state to an excited state. 
When there are many (or even an infinite number of) 
excited states, our considerations can be generalized for 
the analysis of the excitation. In the following we shall 
discuss, then, only the excitation between two quantum 
states. The unperturbed eigenstates of the system are the 
ground state |^)e _lWit and the excited state \ipf)e~ luJft . 
(That is, the eigenstates of the quantum system when 
there is no electromagnetic field due to the blue sheet. 
This could be, for instance, the eigenstates of the sys- 
tem in its original orbit around the BH before the jump 
in, or its eigenstates when the system crosses the event 
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horizon.) Namely, the obvious evolution in time is given 
by the standard oscillatory dependence, and then \ip n ) is 
independent of the time, where n = i, /. The perturbed 
wave function would be given then by 



„ final 



-%UJ n t 



(32) 



where a n are the expansion coefficients. The system is 
taken initially at its ground state Therefore, we 

take a' mtlal = 1 and a™ ltlal = 0. Assuming small transi- 
tions (an assumption which will be justified a posteriori) , 
we take af nal w 1 too. The system is perturbed by the 
pulse of the electromagnetic field. In the Coulomb gauge, 
which is consistent with our treatment, the interaction 
Hamiltonian is Ti = (e/m)A • p, where we keep only lin- 
ear terms in the perturbation, which is assumed to be 
small. Here, A is the electromagnetic vector-potential 
written in the Coulomb gauge, and p is the 3-momentum. 
(The explicit form of A is not important for our needs 
here, although it can be found from E = — <9A/ dr — V0 
and B = V x A, where <f> is the temporal component of 
the four-potential.) It can be easily verified that any vec- 
tor potential of the form A = A(v)e y , where e y is a unit 
vector in the y direction, is consistent with these expres- 
sions and with Eq. (p2|). It can be further shown that 
this form for the vector potential is consistent also with 
the Coulomb gauge condition, namely, with V • A = 0. 
In what follows, we take the temporal component of the 
four-potential to vanish. This can be done consistently 
with the Coulomb gauge. The vector potential can thus 
be written explicitly as A = — J E(r') dr' . 

We are interested in the effects of the blue sheet, 
namely, with the effects of the pulse of the electromag- 
netic field on our system. Therefore, we shall assume, for 
the sake of brevity and simplicity, that the electromag- 
netic pulse ends at the CH, or, in other words, that for 
each point of the system the electric and the magnetic 
fields vanish for positive proper time. More accurately, 
if the system's spatial extension in the z direction is 8z 
from the center of the system, we shall examine the exci- 
tation of the system at proper time (of the system's cen- 
tral point) t > 6z > 0. More conveniently, we shall look 
at the system's state at V > 0, when there is no electro- 
magnetic field associated with the blue sheet. [Of course, 
as physics beyond the CH is as yet unknown, we do not 
suggest here that there are no perturbing fields on the 
other side of the CH. Our point here, is that for the sake 
of the calculation of the blue-sheet effects, the specific 
form of the fields beyond the CH is irrelevant. Moreover, 
our choice here (to set the electromagnetic field equal to 
zero for positive v) is no worse than any other choice 
(in view of the present knowledge of the physics beyond 
the CH).] We shall calculate the energy absorbed by the 
system, as a measure for the strength of its excitation 
due to the divergent electromagnetic field. By first order 
time-dependent perturbation theory a^ nal is given by 



"Km I ^/l A 'Pl^>e' 



dr\ 



(33) 



where f2/i = Wf — u>i. Integrating by parts, we find that 
~ final - ' 1 -e^^/IA-pl^Hr) 
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to h£l fi 
e 1 
to hflfi 
e 1 

TO TlU.fi 

e 1 

to hfifi 

a n + a fi ■ 



(ip f \dA/dt-p\ipi) e in ^dt 
e^(V/|A-p|^)(r) 
(V/|E-p|^)e^*^ 



(34) 



The first term does not contribute to the physical ex- 
citation. The problem with Oj) is that apparently one 
could (incorrectly) infer that even after the perturbation 
vanishes, a^ na[ continues to evolve. To be more spe- 
cific, it looks as though the system is being perturbed 
(and excited) even when there is no perturbation at all. 
This is, of course, an erroneous result. We shall resolve 
this "paradox" in Appendix B, and conclude that afi 
describes a pure gauge distortion of the wave-function, 
which can be set equal to zero by a proper gauge trans- 
formation. That is, when we adjust the gauge such that 
A vanishes for r > 0, the above problematic term sim- 
ply disappears. Therefore, we are thus left only with 
the gauge-independent energy absorption £q(r) which is 
given by 
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a fi 
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[ e l ^*EdtYp 
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(35) 



The physical energy-absorption we are interested in is the 
energy absorption just after the system has fully crossed 
the CH, in the meaning described previously. As we are 
interested here in the effects of the blue sheet, let us 
suppose then that after the CH the electric field vanishes, 
and therefore the quantum system does not undergo any 
further excitation. The integrand in Eq. (|3^) vanishes 
for v > 0, and therefore we may change the variables in 
the integration, and get 



£ q (v = 0) 



1 e z 



\(ip f \e 



Ml fi TO 2 

e m f^E(v)dr ) Pu 
1 e 2 
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(36) 



where p y = p • e y . The integral in Eq. (^) is the same 
integral as in the expression for the total mechanical en- 
ergy of the classical oscillator at the CH, namely, the 
integral in Eq. (Al). Therefore, using Eq. (A27) we get 
that to the leading order in (k-Uq) -1 



£ q (0) = 



1 



(21 

1 
2 



2) 2 nn 

-2(22+2) 

K—Uq 



1 / e \ 2 C" 2 , . , 



(37) 



The matrix element in Eq. (|37j) is regular and does not 
diverge on the CH any more than anywhere else. We 
also see the close correspondence between the classical 
and the quantum systems. As the matrix element in Eq. 
( |37| ) has dimensions of momentum squared, we see that 
the energy gap between the states times the probability 
amplitude for the quantum system indeed has the di- 
mensions of mass. In the classical system the absorbed 
energy is equal to the product of half the reduced mass 
of the system and the square of the internal velocity of 
the system. We see, that the two expressions for the 
absorbed energy in the two systems indeed correspond. 
The advantage of the classical treatment is that it does 
not involve perturbation theory. It suffers, though, from 
the fact that it is a classical model, while actual physical 
matter is intrinsically quantum. Thus, the two models 
contribute to the understanding of each other, and aug- 
ment our understanding of the interaction of the blue 
sheet inside a RN BH with infalling objects. 



VIII. DISCUSSION 



Even though we toy-modeled the matter comprising the 
infalling object in a very simplified way, we believe that 
our models capture the essence of the interaction of ra- 
diation with matter. Therefore, we conclude that the 
classical Maxwell radiation created during the collapse 
of the star will not necessarily destroy objects at the CH. 
We find, that the interaction of the blue sheet with phys- 
ical objects is finite. Moreover, for typical parameters of 
astrophysical supcrmassive BHs this interaction is even 
arbitrarily small. In fact, one can diminish the extent of 
the interaction by simply waiting outside the BH before 
jumping into it. This means that if a spaceship is in or- 
bit around the BH, and an astronaut wishes to fall into it 
(hoping to re-emerge in another universe), he should just 
wait in the spaceship, and postpone the beginning of his 
unusual odyssey. According to our analysis, the longer 
he waits, the smaller the interaction, and the safer the 
voyage. 

We should remember though, that throughout this Pa- 
per we have ignored other inherent radiation sources, 
such as the cosmic background radiation. In addition, 
a more realistic treatment of the interaction of the blue 
sheet with infalling objects will have to consider Qed 
effects, which we have ignored here completely. Such 
Qed effects, and in particular effects caused by pair- 
production, are expected to change the general picture 
portrayed by our analysis considerably. These quantum 
effects might be crucial for a more complete understand- 
ing of the blue sheet and its interaction with infalling 
objects. We showed elsewhere, that these Qed effects 
could be fatal for a human-being observer (due to his 
high vulnerability to 7 rays), but typical physical objects 
of similar or smaller size might survive it. These results 
do not provide support to the idea that no continuation 
of the geometry beyond the CH is physically reasonable 
& 



In this Paper we investigated the following question: 
Are physical objects necessarily burnt up by the blue 
sheet inside a BH? This is a key question for a more com- 
plete understanding of the internal structure of BHs and 
for an understanding of the possibility to fall into a BH 
and re-emerge in another universe. To address this ques- 
tion, we analyzed the interaction of the blue sheet with an 
infalling physical object. We first derived an explicit ex- 
pression for the electromagnetic field as measured in the 
rest frame of the infalling observer, and then we modeled 
the interaction of the blue sheet with the observer in two 
ways: a classical model and a quantum model. In both 
models we calculated a measure for the strength of the 
interaction of the blue sheet with the infalling object. 

We have shown, that the divergence of the energy den- 
sity of the radiation (or, even, the divergence of the inte- 
gral of this energy density over proper time) at the CH 
does not mean that any physical object falling into a BH 
will necessarily be completely burnt up by the radiation. 
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APPENDIX A: EVALUATION OF THE 
INTEGRAL 



In order to evaluate the integral 



/ = / e MT E(t') dr' 

' —OO 



(Al) 



let us divide the integration region into three qualita- 
tively different regions, denoted by a, b and c, respec- 
tively. In region a we assume that — 00 < t <C —lu^ 1 . 
Since the variation of E with r is dominated by 1/r, in 
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region a the electric field changes very slowly compared 
with the exponential in the integrand. Hence, we shall 
approximate the evaluation of the integral by taking the 
electric field outside the integration. We thus obtain that 
I a ps e !WT E(t) i '(ioo) . We now evaluate the error associ- 
ated with this approximation: Integrating by parts (in 
region a), we obtain 



1 



or 
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e luJT E(r) 



r E{r) 



dr' 



1 dE 



ioo dr 

Hence, the error in the evaluation of I a is 
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' , 1 dE 
E(t)--—-E(t) 
iuj dr 



1 



(A2) 



(A3) 



dE 



dr 
(A4) 

As we can take (for region a) E(t) C"/(k_t), we get 
that 



I A/a | 



C 



1 



The relative error is, thus, 



AI a 








1 








Ia 




i 

U1T 




LOT 



(A5) 



(A6) 



We choose now for the proper time r at the boundary 
between regions a and b the following value: 



and therefore we obtain 

AI a 



-1/2 



K- 
< 1. 

U) 



(A7) 



(A8) 



[For typical parameters (see below) k_/lu is of order 
10 -29 .] Hence, we find that the approximation we take 
for region a is valid. 

Region c is defined by — lo^ 1 <C r < 0, i.e., by the 
requirement that the electric field varies very fast com- 
pared with the exponential term, and therefore we can 
take the latter outside the integration. Before we per- 
form this explicitly let us deal with region 6, where most 
of the problems lie. 

Region b is in between regions a and c. The main 
difficulty in the evaluation of the contribution of region 
b to the integral (Al) comes from the neighborhood of 
r « — uj^ 1 : In that region neither of the two approxi- 
mations (the one for region a and the one for region c) 
is valid. It turns out, however, that in that region there 
is a different approximation we can use; the logarithm- 
dependent term in the electric field (E2) does not change 



much in comparison with the 1/r term in region b. It 
turns out, that if we assume that lo^ 1 <C M <C — uq 
(which is a very plausible assumption for physical BHs) , 
region b overlaps with both regions a and c. Thus, using 
the three different approximations (for regions a, b, and 
c), we can calculate the integral (Al) for the entire inter- 
val — oo < t < 0. In fact, the matching between regions a 
and b is done automatically because of the combination 
of two facts: first, in region a the integral follows the 
electric field adiabatically, and does not have a 'memory' 
of its values in former times; second, region b overlaps 
with region a. Therefore, the integral (Al) assumes the 
form 



E 

K- 



1 



dr' 



E 



-(21+3) 



(A9) 



where E = C (]r\k^uu~ 1 \ + ift_u ) X ^' J ' . We justify 
the protraction of the lower limit of the integration in- 
terval from —T to — oo by the vanishing value of this 
expression in the limit of r — > —oo. Region b ex- 
tends up to r = r . We take t such that d = uj\t \, 
where § <C 1 is a dimcnsionless constant, to be evalu- 
ated. (In Ref. (2^] we calculate an optimal value of the 
boundary between regions b and c, i.e., we calculate 
However, this optimal value for -d is of no crucial im- 
portance. Yet, it turns out that this optimal value is 
d « -{21 + 3)(k_u /2)- 1 In \(2l + 3)(k_u /2)^ 1 |.) This 
means that both approximations a and b are valid near 
t — t q . We write the integral in Eq. (A9) as 



1 



dr 



1 



d(u>r). 



(A10) 



[In the approximation for region b we suppose that 
uj- 1 < M < -u : we take the BH mass to be 10 9 M Q , 
where M Q denotes the solar mass. (Even though a typical 
mass for a supermassive galactic BH may be "only" about 
10 7 -10 8 M Q , this does not affect our analysis.) This mass 
is equivalent to a time period of 5 • 10 3 seconds (as 1M Q 
is equivalent to 5 usee). We also take the time of the 
jump into the BH to be of the order of magnitude of a 
typical cosmological time scale, e.g., we take u = — 10 9 
years. This means that k_uo is of the order of —6 • 10 12 . 
If we take oj to be of the order of 10 16 sec^ 1 , which is 
a typical angular frequency for atomic processes, we get 
that In |k_cj _1 | « —40. The typical infall time for a BH 
T sa M, and therefore In |k_T| w 1. This means that the 
variation in the logarithmic dependent term throughout 
the protracted region b is negligible in comparison with 
the magnitude of k-Uq, which justifies the approxima- 
tion we made for region b.] Hence, we get that in region 
b the integral becomes: 



, 6 cos ujt 

h = I d(UJT) + I 

LOT 



smwT 



LOT 



d(wr), (All) 
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or, more conveniently, 

h = {h+il 2 ). 
We now calculate each term separately: 



(A12) 



h=f ^dy' + r^f d y^ l -n + Si(y), (A13) 

J-oo y Jo y 2 

where Si(r) is the sine integral defined by Si(z) = 
Jo ^T^ ^ lUl anc ^ y ~ WT b- As m the overlap region be- 
tween regions b and c we have —lot <C 1, the sine integral 
is much smaller than unity, and we thus obtain that 



!■} W — 7T. 

2 2 



(A14) 



Turn now to I\: 



f v cosy' f 1 cosy' , , , f v cosy' , 



h = i n + i u 



(A15) 



It is straightforward to show that the integral In is 
bounded and of order unity. An explicit calculation 
yields^: 

i n = f 1 £^J/ dy , M 3374 (A16) 



We now write the integral I\2 as 
In -- 



» -1 + cosy' , f 1 , 

— ^ — d y + / 3 d y > 



where K = Ci(l) + A w 0.5772, and where we kept only 
the leading terms. We see from Eq. (A19) that If, con- 
tributes to both the real and the imaginary parts of /. 
However, we now explicitly see the reason for which re- 
gion b can not be protracted all the way to y = 0: the 
logarithmic term in Re(J(,) diverges as y — ► 0. In Ref. 
PH we calculate WT},. It is shown there, that when the 
value of y at the boundary between regions b and c is 
taken for the evaluation of the logarithmic term in Eq. 
(A19), the contribution of h to the real part of / is neg- 
ligible in comparison with the contribution of I c to the 
real part of /. Therefore, we obtain that on the bound- 
ary of regions b and c, I is proportional to (K_u ) _( ' 2i+3 - ) 
due to the proportionality to Eq. Let us now evaluate 
the error in our approximation, and thus show that the 
approximation is valid. The exact value of the integral 
(Al) in region b is 



T 



In \k—t\ + — 



-(21+3) 



= ( hl|K_T b | + -K-UQ 



( s -(21+3) ,r b l 



— e 

T 



(21 + 3) 



x [ ln|/t_r| + — K-Uo 



-(21+4) 



dr. 



Our approximation is 



japprox. _ ( m | K Tb | _|_ - K _u 



-(21+3) 



— e 

T 



dr. 



Therefore, the relative error is 



- 1 + cos y' 



dy' + ln|y| 



(A17) 



The integral on the right-hand side of (A17) can be writ- 
ten as: 



- 1 + cos y' 



dy' = A + 0(y 2 ), 



(A18) 



where A = /,(— 1 + cos y')/y' dy'. Again, it is straight- 
forward to show that the integral A is bounded and of or- 
der unity, and numerical calculation yields^ A w 0.2398 
Therefore, we get that 



Tapprox. 

h 



rcxact 
2 b 



^-approx. 



(2/ + 3) 



In \n-Tb\ + \k^u q 



(21 + 3) 



In \uiTb\ 



In |K_T b | + 2«-Uo 



Taking now the optimal value for the boundary between 
regions b and c (see Ref. P2|), we obtain 



h ^ K + ln\y\ + 0(y 2 ) + i 







AI b 


m (21 + 3) 


In | iK^wo| 


0(y)_ 


(A19) 


h 





< 1. 



(A20) 



^ In fact, In = Ci(l), where Ci(z) is the cosine integral 
defined by C\{z) =j + In z + f* [(—1 + cos t)/t\ dt, where 7 is 
Euler's constant [ p8[ . 

" In fact, A = g° =1 (-l)f +1 - 



Thus, our approximation for region b is valid. 

Region c is defined by the requirement that ?9 < lot < 
0. This means that the electric field varies very fast com- 
pared with the angular frequency of the oscillator. There- 
fore, we may take the e tWT term out of the integration in 



2p(2p)! 
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Eq. (Al). The remaining integral is easily solvable, and 
we get 



1 



1 



K - T ' (ln| K _r'| + i K _M ) 2 ' +3 
An explicit integration yields, then 

1 f 1 



dr'. (A21) 



(2/ + 2) K _ | (in| K _ T | + i K _ Uo ) 
1 ] 



21+2 



which for r = becomes 



/ c (t = 0) 



(2? + 2)k_ 



1 



2i+2 f ' 



-(2Z+2) 



(A22) 



Let us now evaluate the error of our calculation for region 
c. The exact integral for r = is 



1 



In \k—T I H — k_u 



-(21+3) 



, I f 1 X -( 2i + 3 ) 

(1 + iur') - In |k-t'| + -k-Uq 

K-r V 2 



/■T=0 


1 








/■T=0 


+i — 


K- J 





1 

•2' 

1 



-(2i+3) 



-(2Z+3) 



dr' 
dr' 
dr' 
dr'. 



Now, the first term in the right-hand side of the last 
expression is the same as our approximation for I c [see 
Eq. (A21)]. Therefore, we obtain for the relative error of 
our approximation for region c: 



^apprc 



^approx. 



(r = 0) 



^f:-°(ln| K _r'| + !«_no) 



-(21+3) 



dr' 



1 (\ \-(2l+2) 
(2i+2)fc_ \2 K - Uo ) 

^(^«o)- (M+ " 



WT t /" 1 

— K_Mq 



(2; 
-1 



rkrr (l K - u o) 



-(21+2) 



2/ + 3 /l 

— K-Uq 



21 + 2 \2 7 2Z + 2 \2 
For the typical parameters we chose we thus obtain 

jexact ^approx 



^approx. 



(r = 0) < 1. 



(A23) 



For the evaluation of Eq. (A23) we again used the results 
of Ref. 0] . From the evaluation of the error associated 
with our approximation for the contribution of region c 



to 7, it is clear that the most dominant error in I c is imag- 
inary. Therefore, we should also verify, that this imagi- 
nary error is negligible in comparison with the imaginary 
part of the contribution of region b to I. Repeating the 
relative error evaluation, we now obtain, 



jappr. 



(£o/*-)Im(4 aPPrOX - 

(•r=0 /, , n 1 \ -(21+3) 

i Tb (In I s_t I + a/c_u J 



0) 



dr' 



£f(5*-«o) 



< 



U)T b 




2(21 + 3) 


ir/2 




7T 



-(21+3) 



-K-Uq 



< 1. 



Collecting the contributions to I from all three regions, 
we find that the dominant contribution (at r = 0) to 
Re(J) comes from region c. This contribution is propor- 
tional to (k-Uo)~( 21+2 ^ . The dominant contribution (at 
r = 0) to Im(7), however, comes from region b, and is 
proportional to (k-Uq)~( 21+3 \ Hence, we find for the 
required integral 



7(t = 0) = Re [1(0)} + iJm [1(0)} , (A24) 



where 



Re [1(0)} 



Ql / x \-(2l+2) 



-O 



(2Z + 2)k_ V 2 
( K -n )- (2l+3) 



and 



Im[/(0)] = ^Q«_«o 



-(21+3) 



o 



(A25) 



-(21+4) 



(A26) 

When r = we thus obtain, to the leading order in 

(K-«o) - \ 



1(0) « - 



a 



(2l + 2)n_ \2 



1 

— K-Uq 



-(21+2) 



C n fl 

H 1 —K-Un 

k-2 \2 



-(21+3) 



(A27) 



Note, that X(t = 0) is dominated by region b [Eq. (j3l|)], 
and X(t = 0) and £ c (r = 0) are dominated by region c 
[Eqs. © and (§(])], as well as £ q (0) [Eq. (§7 



APPENDIX B: THE QUANTUM FORMALISM 

Consider an electric field E x (v) in flat spacetime, which 
vanishes for v < vq. After v — vq the field E x (v) rises, 
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and vanishes again for v > V±. This electric field can 
be obtained from a potential, whose only non-vanishing 
component is A x (v), which vanishes for v < vq, starts 
changing at v = vq, and assumes a final constant value 
Ao for v > V\ . (We note that any such field is consistent 
with the Maxwell field equations.) We notice that 



A 



E x (v) dv. 



(Bl) 



Let us take a quantum system, which initially is in the 
eigenstate We calculate a/ according to first-order 
time-dependent perturbation theory. In the Coulomb 
gauge, in which our vector potential is written, the in- 
teraction term in the Hamiltonian is — j-^A ■ p. We cal- 
culate af for v > v\. More exactly, we assume that 
the system is centered about z = 0, and that its typi- 
cal spatial extension is 8z. We calculate a/ at the time 
v-2 > Ti + 8z, such that for all points of the system 
v > V\. This means that the system has fully crossed 
the region of non- vanishing electromagnetic field. Hence, 
the electromagnetic field vanishes everywhere at the in- 
stant we calculate a/. We now notice that at that instant 
A x — Aq ^ in general. Therefore, the interaction term 
does not vanish formally, and from first-order perturba- 
tion theory we get 

h f = -~(f\A oPx \*y n ^, (B2) 

which generally is non-vanishing. Integration of (B2) 
yields 



«./■ 



i e 



const - 
nm \ 



A 



(f\Px\i)e 



ifta fiT 



(B3) 



(Ao is constant, and therefore we could take it out of 
the matrix element.) This may lead us to a weird (and 
wrong) conclusion: even though there is no electromag- 
netic field, af apparently continues to evolve. The res- 
olution of this apparent "paradox" is as follows. If we 
write the state of the system for v > v 2 in the standard 
gauge appropriate for that state (namely, A x = instead 
of A x = Ao), we find that a/ is constant. [Initially, since 
we started with E — 0, i.e., with a vanishing perturbing 
field, we chose the standard gauge for the description of 
the system, namely, we chose A x — 0. (This is, of course, 
the gauge in which the system's wave- functions take their 
standard form.) Now, when we come to interpret the fi- 
nal state (in which again E = 0), we again use the gauge 
which is most natural for this situation, namely, we again 
take A x — 0. This requires us to transform our original 



results [Eq. ( 36 )] to the new gauge.] 

Let us show this now: The full wave-function is given 
by 



1&(r) = OifrMe-* 1 * T + E a/(r)|/)e- 
{!/)} 



(B4) 



We take a^(r) w 1, which is consistent with small transi- 
tion amplitudes. Substituting (B3) in (B4) we find (for 
v > v 2 ): 



m = ( K) - E w 



i e Ao 
h m iQfi 



(f\P*\i)\f) 



{!/>} 



(B5) 



The constant term in Eq. (B3) may depend on /, and 
therefore was given in Eq. (B5) an appropriate index. 
The gauge transformation 



A™™ = A° ld - A Q 



(B6) 



of the potential leads to the following gauge transforma- 
tion of the wave-function: 



r"' v ' cxp ( j I An<ir ] < 



old 



(B7) 



Since we work to first order in the perturbation, we may 
use 



r 



1 + — Aox ) i/j 



old 



(B8) 



Substituting Eq. (B5) for V we get that 



1 + -A X 




{!/>} 



x e 



{\f)} 



x e 



(I/)} 

and therefore we obtain 



771 



urn 



(B9) 



after keeping linear terms in the interaction only, and 
using the relation 

E ^{f\PMf)=imx\{). (B10) 



{|/» 



fi 



Hence, we find that the transformation to the natural 
gauge (A — 0) removes the second term in the brackets 
in Eq. (B5). We now get the following conclusion: The 
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correct way to understand the result of the perturbative 
calculation is in the gauge where the vector potential 
vanishes. In this gauge, the second term in Eq. (B3) 
vanishes, and we obtain that 

af w = const, (Bll) 

as should be expected (because the electromagnetic field 
vanishes). 

Finally, we find the specific form for the constant term 

(2) 

of Eqs. (B3) and (Bll). To obtain this, we compare aj/ 
[see Eq. @] to the right-hand side of Eq. (B3) and 
obtain that 

w,»ew = _£. 1 f ^ f \E-p\^}e in ^dt. (B12) 
m nilfi J_oo 

After the perturbation vanishes, this expression is con- 
stant. 
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